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We develop a first-order approximation method for tfie influence of spin on the motion of extended 
spinning test masses in a gravitational field. This approach is illustrated for approximately circular 
equatorial motion in the exterior Kerr spacetime. In this case, the analytic results for the first-order 
' approximation are compared to the numerical integration of the exact system and the limitations 

of the first-order results are pointed out. Furthermore, we employ our analytic results to illustrate 
the gravitomagnetic clock effect for spinning particles. 
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' Imagine the motion of an extended spinning body in the exterior vacuum region of an astronomical source. Let 
■ j^fj-v j^g ^j^g energy-momentum of the extended test mass; then, the motion is governed by the four equations that are 
p\| , given by the dynamical law T^'^. j^ = 0. Consider a representative point inside the extended test body (the "center of 
00 ' mass") such that u^' — dx'^/da is its four- velocity vector and a is the proper time along its worldline. The equations 
of motion can then be expressed relative to the chosen worldline. This has been accomplished in an elegant manner 
by Dixon [1]. Using Synge's world function [2], Dixon has defined the infinite set of multipole moments of T^i, in a 
way that is qualitatively similar to the standard nonrelativistic theory and thereby expressed the equations of motion 
^ as [1] 
Oh 

g: ^ = -2^^c./3<5"'' + ^^ (i) 

^ : ^ = P^u^-P'^u^^+T^r (2) 

.!__( \ Here, the momentum vector P^^ and the spin tensor S^^'^ are the first two moments of T^^^ , while the Dixon force 
I and torque T'^" are given in terms of the quadrupole and higher moments of the test body. It is necessary to add a 
?H supplementary condition to equations (l)-(2) in order to fix an appropriate center-of-mass trajectory. It turns out 
that for an extended body the appropriate condition is [3] 

S^'" P^ = 0. (3) 

The work of Dixon provides a natural generalization of the previous results of Mathisson [4] and Papapetrou [5] that 
were limited to pole-dipole particles. 

In the present paper, we are interested in the dynamics of an extended body of mass m and spin s in the field of 
a rotating mass M m; therefore, we consider equations (l)-(3) without the quadrupole and higher moments, i.e., 
JF^ = and T^'^ — 0. The resulting reduced Mathisson-Papapetrou-Dixon ("MPD") equations then lead to natural 
definitions for the mass m and spin s of the extended test particle that are preserved throughout the motion. We 
assume that the metric tensor has signature +2 and we use units such that c = G = 1, unless otherwise specified. It 

1/2 

then follows from Eqs. (l)-(3) that in this case the mass of the extended particle is given by m = (— P^) ' , which 
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is a constant of the motion. The spin vector is then defined by 

= -^e^upaP'' SP^, S^" = -e^^P'^PpS., (4) 
2m m 

where e^vpa is the Levi-Civita tensor given by Cfj^vpa = \f—g(ypvpa in terms of the alternating symbol with 0-0123 = 1- 
It follows from the reduced MPD equations that = S^^S^ = \ S^^ S^" is a constant of the motion. Moreover, 
the requirement that the dipole force be much smaller than the; monopole force implies that s <C ?7i r, where r is 
the distance between m and M; that is, the M0ller radius of the extended test particle (s/m) must be very small 
compared to its distance from the source. 

In a recent paper [6], an approximation scheme has been introduced to deal with this extended pole-dipole system 
in most astrophysical situations. This scheme is based on the circumstance that — is small and of order 

(M/r) [s/{mr)]^ <i; 1. In this approach, it then turns out that to first order in s/{mr) one can assume that P'^ and 
are parallel, namely, P^" « mu^; therefore, it follows from Eqs. (l)-(3) in this case that 

D7/'* 1 „ nsf^"^ 

- -^R^^^^<S'^^. ^«0, 5,.<.0. (5) 

The second section of Ref. [6] should be consulted for the details of the derivation of these equations. 

In the absence of spin, system (5) simply reduces to the geodesic equation for the motion of the center of mass of 
the spinless extended particle, as expected [7]. In Sec. II, we solve system (5) assuming that in the zeroth order of 
approximation the worldline for a spinless particle is a timelike geodesic. The motion away from the geodesic is then 
of first order in s/{mr) in our approach. This method is employed to study approximately circular motion near the 
equatorial plane of a Kerr source of mass M and angular momentum J = Ma in Sec. III. That is, we assume that 
in the absence of spin, the test mass follows a stable equatorial circular geodesic orbit in the exterior Kerr spacetime. 
Our approximate analytic results based on system (5) are compared with the numerical solution of the extended 
pole-dipole system, Eqs. (1) (3), in Kerr spacetime in Sec. IV. Section V is devoted to the gravitomagnetic clock 
effect. A brief discussion of our results is contained in Sec. VI. The appendices contain further developments of our 
main results. 



II. FIRST-ORDER APPROXIMATION 



Consider the motion of an extended spinning test particle in a gravitational field. Let w,^ be the four-velocity of the 
particle in the absence of spin and A^a be the orthonormal tetrad frame that is parallel propagated along the reference 
geodesic a;^(r), where u^" = dx^/dr and t is the proper time of the spinless particle. More specifically, we imagine 
that along the unperturbed reference geodesic worldline x^^{t), the Mathisson-Papapetrou spin-curvature interaction 
is "turned on" at r = and the spinning particle then follows .x^(t) for r > 0. Let Sx'-^{t) = .x^(r) — .x^(r) be the 
deviation between the neighboring worldlines at the same time r > such that {Sx'^) = 0. We will work to linear 
order in the small quantity s/ (mr) <^ 1. Here Sx^^{t) is a vector field defined orthogonally along x^{t) such that it 
connects the reference geodesic to the path of the spinning particle in a way that amounts to a unique identification 
of events along the perturbed trajectory using the proper time of the geodesic worldline x^^. It follows that we can 
write 

fe'^ = X\t)X'^u (6) 

where X(t) can be determined by the equations of motion (5). To this end, we establish a Fermi normal coordinate 
system along x^^{t) based on the local frame X^^& and solve the equations of motion in this Fermi coordinate system. 
The procedure that we follow is described in detail in Appendix A. Here we provide a general summary of our 

approximation method. 

To linear order in our approximation scheme, it follows from Eq. (5) that the spin tensor S^^"^ is parallel propagated 
along the unperturbed trajectory and S^i, = 0. One can then write 

S^"' = A^jA'^jS'^ (7) 

where S^^ are constants of the motion. Moreover, it is simple to show that in the linear approximation scheme under 
consideration in this paper, S^u^ = 0, so that for Sq. = S^X^Q: we have -Sq = and 



Si — — £ijk S^^ . 



(8) 



3 



Here, we have used the fact that e^^pcr A^qA'^ jA'^^A'^g = 1, since egjgg = 1; hence, det (A'^a) = {—g)~^^'^, which is a 
consequence of the orthonormahty of the tetrad frame. 

Consider next the components of the spacetime curvature tensor measured by a free observer moving along the 
reference geodesic. These components, iJ^^pg^ A^q,A''^A''^^A'^^-, can be expressed as a 6 x 6 matrix (TZai:), where the 
indices A and S range over the set {01, 02, 03, 23, 31, 12}. One can then write 

^=(7?t2)' (9) 

where £, 7i and V arc 3x3 matrices representing respectively the electric, magnetic and spatial components of the 
Riemannian curvature. In general, £ and T) are symmetric and Ti is traceless. In a background Ricci-flat spacetime, 
V = —£, £ is traceless and H is symmetric {H = H^) . In this case, the curvature is completely determined by its 
electric and magnetic components. These are given in general by 

£jj = R^vpa A'^g A'^j A''q A'^j, (10) 

and, as described below, play important roles in our approximation scheme. 

It is shown in Appendix A that the worldline of the spinning particle is given, in our linear approximation, by 



<(t) = xnr)+XXT)A^-, (12) 



where X{t) is a solution of the system 



V\ (13) 



dV^ 1 

— +%(t)X^ = -n-.s{T)S\ (14) 

with the boundary conditions that at r = 0, X{0) = and V(0) = 0. The unique solution of this system is thus 
due to the spin-curvature force folded together with the tidal influence of the background gravitational field. In fact, 
system (13) and (14) can be solved with the method of variation of parameters [8] once the fundamental solution of 
the homogeneous part is available. In the absence of spin, the system simply reduces to the Jacobi equation in the 
Fermi coordinate system; therefore, system (13) and (14) can be solved once the Jacobi fields along the path of the 
reference geodesic are completely known. The general solution of system (13) and (14) is discussed in Appendix A. 
In the following section, we apply this method to Kerr spacetime; that is, we solve Eqs. (13) and (14) for the case of a 
spinning test particle that follows an approximately circular equatorial orbit in the exterior Kerr field. The Kerr field, 
which represents the exact gravitational field of a rotating black hole, has been chosen here for the sake of simplicity; 
in principle, one could employ any other exact solution of the gravitational field equations describing the exterior field 
of a rotating mass [9]. 



III. KERR SPACETIME 



Consider circular geodesic orbits about a Kerr source of mass M and angular momentum J. In standard Boyer- 
Lindquist co-ordinates {t,r,0,(p), the geodesic equation for a circular equatorial orbit of fixed "radius" r > 2M and 
6 = Tr/2 reduces to 

t = -^{l+anK)T, <i> = j^^KT, (15) 
where a = J/M > and we have chosen boundary conditions such that t = (^ = 0atT = 0, 



and 



N = ^ll-^ + 2anK (16) 



M 
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is the Kcplcrian frequency of the orbit. Here the upper sign indicates an orbit that rotates in the same sense as 
the source, while the lower sign indicates a counter-rotating orbit. The energy E and orbital angular momentum L 
associated with these circular orbits are given by 

+ (IS, 

A detailed discussion of Kerr geodesies can be foimd in Refs. [10, 11] and references therein. Furthermore, the 
motion of spinning test particles in Kerr and other astrophysically interesting spacetimes has been studied by a 
number of authors using other approaches, as can be seen from a perusal of Refs. [12-23]. 

It is possible to set up an orthonormal parallel- propagated frame A^^ along the geodesic worldline of the spinless 
test particle such that A'^q = dx'^/dT is the particle's four-velocity vector and A''?, i = 1,2, 3, are unit gyro axes that 
form the particle's local spatial frame. It follows that in {t, r, 9, (j)) co-ordinates [24] 

= ^- A sine, A cose, 0, sine^ , (20) 

A'^2 = (o,0,i,0^ , (21) 

A'^3 = cose, A sine, 0, ^ cos 9^ , (22) 

where e = Okt and A is given by 



^=\/l- — + 72- (23) 
Furthermore, the electric and magnetic components of the spacetime curvature are given by [24] 














h 






k2 




H = 


I 





I 













h' 





£ = I k2 ] , n = \ h h' \ , (24) 
\k' ks J \ 

where k2 = — (fci + fcs) is a constant given by ^2 = (37^ — 2) and 

ki = A; (1 - 37^ cos^ e) , fca = A; (l - 37^ sin^ 9) , k' = -372/0 sine cose, (25) 
h = -37^/3 A; cose, h' = -37^/3 A; sine. (26) 



Here 



^ = ^ = ^K, (27) 



while /? and 7 constitute a Lorentz pair, 7 = — (3^, given by 



The results presented here for the electric components of the curvature tensor are consistent with the work of Marck 
[25]. 

We can now employ these results in order to find the perturbed orbit x^{t) = {ts,rs,dsj'f>s) given by Eq. (12), 
namely, 

^ ^ + sine + Z cose), (29) 

N rA^ ' ^ ' 

= r + A{X cos® + Z sin e) , (30) 

TT Y 

Os = ^ + -, (31) 

= ^^r-h^ (-Xsine + Zcos9), (32) 
N rA 
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where X = (X, Y, Z) and the tetrad frame given by Eqs. (19)-(22) has been used. It is useful to write in terms of 
spherical polar co-ordinates with respect to the local tetrad frame as 



Si ~ s sind cos 



^2 = 



COST?, 



= s sin?? sin ip, 



(33) 



so that s cosd is the component of the spin vector along the rotation axis of the Kerr source. Then system (13) and 
(14) takes the form 



d^Y 
d^Z 



(37^ cos^ e - 1) X + 37^ sin 9 cos 9 Z + 37^ /? cos -dj cos 9, 

(2 - 37^) y - 37^ /? sint?) cos {& - if) , 

37^ sin 9 cos e X + (37^ sin^ 9 - l) Z + 37^ /3 (-^ cos sin 9, 



(34) 
(35) 

(36) 



with the boundary conditions that X = and X = at 9 = 0. Here an overdot denotes differentiation with respect 
to 9. Let us now consider a rotation about the vertical axis by an angle 9 so that the rotated coordinate axes 
correspond to the radial, vertical and tangential directions; that is, Y' = Y and 



X' = Xcos9 + Zsin9, Z' = -X sin9 + Z cos9. 



(37) 



We note that a;^ can be conveniently expressed in terms of the rotated coordinates {X' ,Y' , Z'), as they appear 
explicitly in Eqs. (29)-(32) for the perturbed orbit. 

It is possible to show that under rotation (37), Eqs. (34)-(36) take the form 

X' -2Z' -3"/^ X' = :i-f^ (3(— cosA , (38) 

Y' + (37^ - 2) r' = -372 p sin^) cos {Q - <p) , (39) 

Z' + 2X' = 0, (40) 
with the boundary conditions that at 9 = 0, X' = and -X" =0. The unique solution of this system is given by 



X' = — 7^/3 fi- cos i?) [1- cos (p 9)1, 
Y' = \ sini?) 



cos i-p cos (C 9) + - sin ^p sin (^ 9) — cos (9 — ^p) 



Z' = -^j^p(—cosA[pe-sm(pe)]. 



Here 



\/4-372, C = V3? 



(41) 
(42) 
(43) 

(44) 



are such that pIOrI and C|^k| are the proper radial and vertical epicyclic frequencies, respectively. 

It follows from these results that the orbit of the spinning particle in the standard Kerr coordinate system is given 

by 



ts = 



1 + a Ok 
N 



T + 



■r 



/3(^cosi?)^[sin(p9)-p9], 



TT 1 

2 [3 \mr 
Ok 6 2 



sin?? 



cos (p cos (C 9) + - sin ^p sin [Q 9) — cos (9 — <p) 

E 



^(-^r cos^)^[sin(p9)-p9] 



(45) 
(46) 
(47) 
(48) 



The perturbation in (ps is simply E / L times the perturbation in ts- In Eqs. (45) (48), the temporal parameter r can 
be replaced by a, which is the proper time of the perturbed worldline, since r = cr + 0{s^) by Eqs. (A.7)-(A.9) of 
Appendix A. It is interesting to note here the role of spin length scales 



s* = {—] COSI?, s = ( — ) sini? 
Vm/ Vm/ 



(49) 
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in the worldline of the particle; the perturbations in tg, rg, and (j)g are proportional to s*, while the polar motion away 
from the equatorial plane is proportional to s. 

Equations (45)-(48), involving the first-order perturbation terms proportional to s/{mr), should be viewed as 
containing the zeroth and first-order terms of a perturbation expansion in powers of a sufficiently small parameter 
s/{mr) <C 1. Moreover, for a black hole (a < M), (3 in Eqs. (45)-(48) is positive for co-rotating orbits and negative 
for counter-rotating orbits, while 7^ monotonically decreases from 4/3 at the last stable circular orbit to unity at 
infinity. Similarly, |/3| decreases monotonically from 1/2 to zero over the same range; in fact, for r ^ M and r ^ a, 

It is interesting to mention here the behavior of /? for a > M. As before, f3 is negative for counter-rotating orbits, 
while for co-rotating orbits can be positive or negative. In fact, for a co-rotating orbit in this case (a > M), /3 is 
positive for r > a? /M, negative for r < a? /M, and zero for r = a^/M. 

In Eq. (46), it is important to note the periodic motion in the radius: it oscillates between r and r + 6 ^4 7^ /3 s* / 
with a proper period of 27r/ (p|riK|)- For the motion of the Earth about the Sun, r ^ M and the M0ller radius of the 
Earth is about 200 cm, so that the amplitude of this periodic variation in the astronomical unit amounts to 

— « - I^kIs* ~ 10-", (50) 
re 

which is a few orders of magnitude too small to be measurable at present. In this connection, it is interesting to note 
that a secular increase of the astronomical unit, amounting to about 10^ cm per century, has been recently reported 
based on the analysis of radiometric data [26, 27]. Furthermore, the polar motion about the equatorial plane as 
given by Eq. (47) consists of a superposition of two frequencies: jf^Kl and C|Ok|. For r ^ M > a, this motion in 
the polar direction exhibits a beat phenomenon. That is, the "fast" harmonic oscillation at essentially the Keplerian 
frequency |0k| is modulated by a "slow" beat frequency « i^g/^, where iVg = 3GM|f2K| / (2c^r) is the "geodetic" 
(i.e., de Sitter- Fokker) precession frequency for the unperturbed orbit. The net amplitude (in radians) of the polar 
motion for r » M > a is independent of c and is given by s/VGMr. Nevertheless, this maximum amplitude is 
a relativistic effect and builds up over a "long" timescale; this circumstance is reminiscent of another phenomenon 
related to relativistic nutation [28]. For the Earth, the yearly spin-induced up and down polar motion about the 
ecliptic has a net amplitude of about 10^ cm and beat period of 10* years; that is, the amplitude of the polar motion 
away from the ecliptic develops gradually and reaches a maximum of about 10 kilometers over a period of about 25 
million years. 

It would be very interesting to search for these radial and polar motions in the strong-field regime close to black 
holes or in relativistic binary systems. 

Let us now compare and contrast these analytic (but approximate) results with the numerical analysis of the exact 
equations in Sec. IV. 



IV. NUMERICAL INTEGRATION OF MPD EQUATIONS 

To determine the domain of validity of our linear approximation scheme, it is interesting to make a numerical 
comparison of our results with the full MPD equations, which must be written in a form that is suitable for numerical 
integration. A direct comparison is possible once the reduced MPD equations are expressed in terms of proper time 
a as [29] 



da 



dS°^ 
da 



+ ^R"(3p.e^\.S''P''ut (51) 



F fll^ Ug 



Ryffpa eP\, ) P° 



(52) 



da - - V + 2 + i R,.,. S^^ Sp^ ) ' ^^^^ 

where M = —u'^Pa/w? is a normalization parameter such that g^i^u'^u'^ = —1 throughout the spinning particle's 
motion. 
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We wish to integrate these equations numerically in Kerr spacetime with the same boundary conditions as in the 
linear approximation. Thus, we assume that at cr = the particle starts from x'^{a = 0) = (0, r, 7r/2,0) in Boyer- 
Lindquist co-ordinates. Moreover, we choose initial conditions for the linear momentum such that P'^{a = 0) = 
mu'^{a = 0), where 

<(C7 = 0) = A% (54) 

is given by Eq. (19) and corresponds to a stable circular equatorial orbit around the Kerr source. The initial conditions 
for the components of the spin vector are given by 

S>'{a = {)) = X^'&{(j = 0)S". (55) 

Using Eqs. (20)-(22) and (33) with Sq = 0, we find 

^/// f mLs . . _ ms* mEs . \ 

S^{a = 0) = I — — smip, mAs cosip, , — — sm(p I (56) 

in Boyer-Lindquist co-ordinates. In this way, Eqs. (45)-(48) should provide the linear approximation to the full 
solution of Eqs. (51) (53), since, as shown in Appendix A, r = cr + O(s^). The results of the integration of the 
reduced MPD equations for the specific case of a prograde orbit with m = 10~^M, iD = ip = n/A, a = 0.5M, and 
r = lOM are presented in Figures 1-3. Specifically, the spherical polar co-ordinates (r, 6, (f)) characterizing the orbit 
of the spinning particle in the standard Boyer-Lindquist co-ordinates are plotted versus proper time a in Figs. 1-3, 
respectively. 

The results of the numerical integration of the reduced MPD equations confirm that Eqs. (45)-(48) provide the 
proper linear approximation to these equations. It is important to remark that these linear terms primarily exhibit the 
results of spin-orbit coupling (cf. Appendix B); that is, the rotation of the central source, though possibly important, 
is not crucial for the existence of the main phenomena associated with the motion of the particle in the first-order 
approximation. However, the rotation of the source could play a significant role in the difference between prograde 
and retrograde motion as in the clock effect discussed in Sec. V. 

Beyond the linear approximation, it appears from Figs. 1 3 that for radial motion nonlinear spin terms introduce 
a complex beat phenomenon involving at least two periods. Of these, the higher frequency is about the same as the 
beat frequency in the case of polar motion. Furthermore, nonlinear spin effects are responsible for the reduction of the 
beat amplitude in the polar motion. On the other hand, the azimuthal motion appears to be essentially unaffected 
by the nonlinearities; therefore, only the graphs for the highest magnitude of spin considered here are presented in 
Figure 3. 



V. GRAVITOMAGNETIC CLOCK EFFECT 



It follows from the results of the previous section that the analytic formulas involving s j (mr) to first order are a 
good approximation to the exact solution of the reduced MPD equations: to go further, higher-order terms in spin 
must be taken into account. On the other hand, the analytic formulas can be safely applied within their theoretical 
limits of validity to illustrate a general feature of motion about a rotating source, namely, the gravitomagnetic clock 
effect [30 33]. In its simplest form, the effect involves a spinless "clock" on a circular orbit in the equatorial plane 
of a rotating source. The difference in the proper periods of two clocks on the same orbit but moving in opposite 
directions is given by t+ — r_ « Air J/{Mc^) for r ^ 2GM/(?. Here r+ (t_) is the proper period for the prograde 
(retrograde) motion. To lowest order, this remarkable effect is independent of Newton's constant of gravitation G and 
the radius of the orbit. Moreover, the prograde motion is slower than the retrograde motion. For satellites in circular 
orbit about the Earth, r+ — t_ « lO"*" s; however, measuring this effect will not be simple [31]. 

The purpose of this section is to show that for a spinning "clock", the analogous result to first order in s/{mr) is 

based on azimuthal closure involving Eq. (48). This constitutes a generalization of previous results [18, 19, 34] that 
involved only s ■ J = ±1. 

Imagine a spinning particle revolving around a Kerr source in accordance with Eqs. (45)-(48). After an integral 

number p of complete orbits with p ^ 1, the change in (j)s over a period of time pr^ (pt-) for the co- rotating (counter- 
rotating) clock is 2-^ p {—2TT p) . For p^ 1, we can drop the sin(/!j8) term in Eq. (48) in comparison with pQ; then. 



(57) 
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FIG. 1; The orbital radial co-ordinate rs(a) for various choices of s/{mr), where m = 10 ^ M, r = 10 Af, ■& — ip — 7r/4, and 
a = 0.50. 



we find that the resulting expressions for and t_ are such that the one for changes over to the one for t_ if 
a —> —a and s —s. To first order in a/r and s/{mr), we have 



T± = CTk± — la + 6s j , 



where 



B 



1-^, C 
r 



3M 



D = A 1 



6M 



(58) 



(59) 
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FIG. 2: The orbital polar angle 6s{o-) for various choices of s/{mr), where m — W ^ M, r = 10 M, ■§ — ip — 7r/4, and a = 0.50. 

Tk — 27r/ 117x1 is the Keplerian period and s* given by Eq. (49) is proportional to the component of the specific spin 
of the "clock" along the axis of rotation of the source. Let us note that t+ — t_ computed from Eq. (58) for r ^ 2M 
would result in Eq. (57). 

It is interesting to remark here that in the equatorial plane of Kerr spacetime, 



t+ — i_ = An a, 



(60) 



which is another manifestation of the special gravitomagnetic temporal structure around the rotating source. Here 
i+ and t- are the periods of geodesic motion around the same circular orbit according to the static inertial observers 
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(a) s/(mr) = IQ-^ (MPD) 



(b) s/{mr) = 10 ^ (Linear) 



FIG. 3: The orbital azimuthal angle cj}s{cr), where m = 10~^ M, r = WM, ■& = ip = tt/4, and a = 0.50. 



at spatial infinity. Substituting Eq. (58) in Eq. (45), we find along the same lines as before that the corresponding 
result for the motion of a spinning particle to first order in s/{mr) and a/r is t± = Tk ± Svr {a + 6s*C^/Z3^), so that 

t+-t- = 4^ + 6s* . (61) 

These results are in qualitative agreement with previous work in this direction [18, 19, 34] that has been restricted 
to circular orbits for the motion of the spinning particle around the Kerr source with the corresponding limitation 
that the spin of the particle has to be parallel or antiparallel to the Kerr axis [15]. In this connection, it is important 
to remark here that Eqs. (45)-(48) explicitly forbid a circular equatorial orbit with cosd = ±1. Thus, a direct 
quantitative comparison with previous work is not possible in this case. The contribution of spin to the clock effect 
is further elucidated in Appendix B. 



VI. CONCLUSION 



We have developed in this paper a general first-order approximation scheme for the effect of spin on the motion of 
an extended spinning test particle in a gravitational field in accordance with the MPD equations that are based on 
the multipole expansion method. We have neglected the influence of quadrupole and higher moments on the motion. 
For the astrophysically interesting case of motion in the field of a central source, we must have s/{mcr) <^ 1, where 
s/{mc) is the M0ller radius of the extended test particle. Thus the influence of the spin of the particle on its motion is 
treated to linear order throughout. It turns out that in this approximation the spin is parallel transported along the 
nongeodesic path of the test particle. To illustrate our general approach, the motion of spinning test particles along 
nearly circular equatorial orbits in the exterior gravitational field of a Kerr source has been investigated in detail and 
the results of the first-order approximation scheme have been numerically compared with the corresponding solution 
of the full MPD equations. As expected, the results are essentially identical for s/{mcr) sufficiently small compared 
to unity; moreover, we have described the deviations that appear when s/{mcr) is not so small. 

The astrophysical implications of our results have been briefly mentioned in Sec. III. For the motion of the Earth 
about the Sun, for instance, the spin-induced deviations from Keplerian motion are too small to be observationally 
signiflcant at present; however, the spin-dependent terms could become important in binary pulsars [35, 36]. Fur- 
thermore, we have elucidated the gravitomagnetic clock effect for spinning particles in accordance with our flrst-order 
approximation method. 
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APPENDIX A: MPD EQUATIONS IN FERMI NORMAL COORDINATES 



Consider a Fermi normal coordinate system that is constructed about the reference geodesic a;^(r) and is based 
upon the paraUcl-propagatcd orthonormal frame A^^. Specifically, along the worldline of the spinning particle a;^, let 
us imagine an event Qs with Fermi coordinates = [T, X) ; then, there exists a unique spacelike geodesic of proper 
length I that connects Qs orthogonally to the reference geodesic at the event Qo with proper time r such that 

T = T, Xi = In^Xf^i, (A.l) 

where = {dx^ /dl)Q^ is the unit vector tangent to the unique spacelike geodesic at Qo such that n^w'^ = 0. Thus 
the spinless reference particle following a;^(r) permanently occupies the spatial origin of this Fermi coordinate system 
and has Fermi coordinates (T, 0), where T = r. The spacetime metric in Fermi coordinates is given by 

500 = -1 - i?6?o,-(^) X'X^ + ---, (A.2) 
90i = -^%fc(T)X^X^ + ... , (A.3) 

9ij = Sij - I R,k^({T) X''X' + ---, (A.4) 

where -R^^^^, given by the projection of the Riemann curvature tensor on the tetrad frame of the spinless reference 
particle, is in fact the Riemann tensor in the Fermi coordinate system evaluated at its spatial origin. The Fermi 
coordinate system is admissible within a cylindrical spacetime region around the worldline of the reference geodesic 
such that 1^1 < TZa{T), where TZa is a certain minimum radius of curvature of spacetime. 

The Mathisson-Papapetrou-Dixon equations of motion of the spinning particle within the framework of our first- 
order approximation scheme (based on s <C mr) can be expressed in the Fermi coordinate system as 

= = -^R^^^U-rS'^^ (A.5) 

^— = 0, S^^U'' = 0, (A.6) 
where Lf^ = dX^^/da is the four-velocity of the spinning particle. Let us write this as 

W^TiLV). r = f. F=^. (A^7) 

Then, 

= ~goo-2go^V' -g,,V'V^ > 0, (A.8) 

since f/'^ is a timelike unit vector. To characterize the order of various quantities in accordance with our perturbation 
method, it is convenient to write 

X{T) = 0{s), V{T) = 0{s), ^ = l + 0{s'), (A.9) 

and so on; in particular, it follows that a = T + O(s^). We have already discussed the consequences of Eq. (A.6) in 
Sec. II; therefore, we concentrate here on Eq. (A.5). We find [6, 37] 



where 



A' = 0, A' = ^i?Q^,^(T)5^-^ (A.11) 



We note that the Mathisson-Papapetrou acceleration is given by 

1 



Hi-i (A.12) 
m 
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based on the definitions introduced in Sec. II. Moreover, of the connection coefficients in Eq. (A. 10) only Fqq = £ij 
makes a non-negligible contribution. Thus taking due account of our first-order approximation scheme, we find that 
Eqs. (A.10)-(A.12) reduce to 

^+e,iX^ = (A.13) 

with the boundary conditions that X = and V = at T = 0. We now turn to the general solution of Eq. (A.13), 
which may be expressed as 



where / is the 3x3 unit matrix and 

Let V be a general solution of the homogeneous (Jacobi) system with x = 0; then. 



*=fjV X=f°. ). (A.15) 



V' = l^Ci^i, (A.16) 



where the Cj, i = l,2,---,6, are arbitrary constants that physically correspond to the initial position and velocity of 

a free particle following the Jacobi equation and V-'i, i = 1, 2, • • • , 6, form a fundamental set of solutions of the Jacobi 
system. According to the method of variation of parameters [8], we seek a solution of the inhomogeneous Eq. (A. 14) 
by assuming that Cj Ci{T), i.e., we let 

6 

* = ^a(T)Vi. (A.17) 

1=1 

Substitution of Eq. (A.17) in Eq. (A.14) results in 

E^^^ = X. (A.18) 



i=l 



Consider a 6 x 6 matrix $ that is the fundamental matrix solution of the homogeneous Jacobi system and contains 
V'l) V'2j • ■ ■ > V'e as its column vectors. Then = $C, where C is a column vector with Ci, C2, • • • ,Cq as its elements, 
and Eq. (A.18) can be written as 

§ = ^-'X. (A.19) 
With the initial conditions that ^'=C = OatT = 0, the solution of Eq. (A.14) is 

*(T) = $(T) / ^-\T')x{T')dT'. (A.20) 
Jo 

In principle, the trajectory of the spinning particle can thus be determined in the Fermi coordinate system. To 
express this trajectory in terms of the original background coordinate system, it is in general necessary to have the 
explicit coordinate transformation between the two systems of coordinates. This turns out to be possible only in very 
special situations [38, 39]. However, our first-order approximation scheme makes it possible to proceed as follows. 
The deviation vector Sx'^ may be written as Sx'^ Ri In^, since I <C TZa- Therefore, regarding Sx^ as a vector field 
along the reference geodesic, we can write Sx'^ = X'^{t) X^i, where the Fermi temporal coordinate in the solution of 
Eq. (A.13), namely, X{T), has been replaced by r. It follows that a;^ = a;'* + X' A''?; in this way, we recover Eqs. (6) 
and (12) of Sec. II. 
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APPENDIX B: HEURISTIC INTERPRETATION OF THE CLOCK EFFECT FOR SPINNING 

PARTICLES 



The purpose of this appendix is to provide an intuitive physical understanding for the appearance of terms pro- 
portional to s* in the discussion of the gravitomagnetic clock effect in Sec. V. The existence of such a term was first 
pointed out in Ref. [18] and has been further studied in [19] and [34]. The approach adopted here is based on the 
gravitoelectromagnetic ( "GEM" ) analogy; we follow here the basic conventions of Ref. [40] . 

Consider a spinning particle in orbit about a central nonrotating source. In the rest frame of the particle, the 
central source revolves around the particle and this motion generates a gravitomagnetic field. The resulting spin- 
gravitomagnetic field interaction, just as in the electromagnetically analogous case of the hydrogen atom, must be 
corrected by taking due account of Thomas precession, which in the gravitational case becomes the de Sitter-Fokker 
("geodetic") precession. The combined interaction due to these terms is of the spin-orbit coupling form [36]. Specifi- 
cally, in the rest frame of the spinning particle, the gravitomagnetic field Bg is proportional to — {v/c) x Eg, where 
Eg = GMr/r^. The corresponding contribution to the Hamiltonian, i.e. the analogue of — /x • B, would be [40] 

-s-Bg = ^J-^^\s-L, (B.l) 

where L = mr x v and is a numerical factor that is expected to be of order unity. Furthermore, the geodetic 
precession frequency is given by 



3 / GM 



2 \ mc^ 

hence the corresponding contribution to the Hamiltonian would be 



(B.2) 



Therefore, the net Hamiltonian for the motion of the particle may be taken to be of the form 

2m r ymc^r'^ ) 

in this simple GEM model. Here, 

^ = 6 + 1 (B.5) 

In Eq. (B.4), the spin-orbit term is a small relativistic correction to the Newtonian dynamics and turns out to be the 
source of the spin-dependence of the clock effect [34] . 

It follows from Hamilton's equations of motion for Eq. (B.4) that 

dr V .( GM 



. - .a ]sxr, (B.6) 

dt m \ mcr J 

dp fGMm\ ^JGM\, / GM \ 

Thus the canonical momentum can be expressed as 

p = mv + ^\^-^jrxs (B.8) 

and the equation of motion of the spinning particle is then given by the substitution of Eq. (B.8) in Eq. (B.7). To 
first order in the spin-dependent contribution, the result is 

fGM\ F , , 

where F is the spin-dependent force given by 

jGM\{Zdr, ^ X 3 , , 1 

^ = ^(^j|-^(^x«)-2(^xs) + -^[(rxt;).5]r|. (B.IO) 
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Consider an initially circular Keplerian orbit of radius r in the equatorial {x, y) plane. At t = 0, the orbit is 
perturbed by the spin-dependent force 

F = ^Ok(^|^) [(scos^)r + 2(sT)z], (B.ll) 

where s = (s, "d, if) in spherical polar co-ordinates and 

r = r[{cosnKt)x + {sinnKt)y]. (B.12) 
Let us write the equation of the perturbed orbit in cylindrical co-ordinates as 

rs = r{l + f), (j>s = ^Kt + q, = rh, (B.13) 
where /, q, and h denote spin-dependent perturbations. Imposing the boundary conditions that at t = 0, 

we find, using the method developed in Ref. [41], that 



= r + ^^^s*{l-cosnKt), (B.15) 
= flKt--^flK s* {flKt - sin ClKt) , (B. 16) 



^ — s [sin {flKt — (p) f^K^ + sin ip sinf2K^] • (B.IT) 
r 

It is interesting to compare Eqs. (B.15)-(B.17) with Eqs. (45)-(48); in fact, to lowest relativistic order, we find from 
the latter equations the same results as Eqs. (B.15) and (B.16) with ^ = 3. Moreover, with (^Q « I^k^, C ~ 1 + (3/2) 
and w GM/{c^ r), Eq. (47) implies that 



TT 3 £ [gin ^f^xi — tp) f^K^ + sin sin f^K^] • (B.18) 



2 2 V c2 r r 

We note that Zg = Vg cos'ds, so that one obtains from Eqs. (45)-(48) that 



3 / GM 

Zs i=a 2 V ~^ — s [sin (Ok^ — ¥') ^^K* + sin sin J^xi] , (B.19) 

which is smaller by a numerical factor of 1/2 than the expression given by Eq. (B.17) with ^ = 3. Nevertheless, it is 
rather remarkable that our simple model, based on the GEM analogy, predicts the main qualitative features of the 
perturbed orbit. In particular, it follows from Eq. (B.16) that 

t± = T^i 477^5*; (B.20) 

hence, for ^ = 3 the spin part of the clock effect is recovered. 

We have thus far ignored the proper rotation of the central source. The rotation of the source would generate a 
gravitomagnetic field and the spin of the test mass naturally couples to this gravitomagnetic field as described in 
detail in Ref. [6]. It turns out, however, that the corresponding dominant "hyperfine" coupling term is independent 
of the sense of the orbit and hence does not contribute to the main spin-dependent gravitomagnetic clock effect under 
consideration here. 

Finally, a remark is in order here regarding the fact that with ^ = 3 and ^ = + 3/2, we have = 3/2, so that 
the s ■ Bg/c and geodetic terms, given respectively by Eqs. (B.l) and (B.3), contribute equally to the net spin-orbit 
coupling term in the classical Hamiltonian (B.4). As noted in Ref. [34], a similar result has been obtained in the 
treatment of a Dirac particle in the Schwarzschild field [42-45] . 
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